Chaotic enhancement in microwave ionization of Rydberg atoms 
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The microwave ionization of internally chaotic Rydberg atoms is studied analytically and nu- 
merically. The internal chaos is induced by magnetic or static electric fields. This leads to a chaotic 
enhancement of microwave excitation. The dynamical localization theory gives a detailed descrip- 
tion of the excitation process even in a regime where up to few thousands photons are required to 
ionize one atom. Possible laboratory experiments are also discussed. 
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I. INTRODUCTION 

The pioneering experiment of Bayfield and Koch per- 
formed in 1974 attracted a great interest to ionization 
of highly excited hydrogen and Rydberg atoms in a mi- 
crowave field • The main reason of this interest is due 
to the fact that such ionization requires the absorption 
of a large number of photons (about 20 — 70) and can be 
explained only as a result of the appearence of dynamical 
chaos and diffusive energy excitation in the correspond- 
ing classical system. Indeed the critical border Cc for the 
microwave field intensity above which classical chaotic 
motion takes place is given by j|] : 



eo > Cc 



1 



(1) 



Here, e and lo are the microwave field strength and fre- 
quency, eo = en-Q and luq — uin^ are the rescaled values, 
and Eq = —l/2nQ is the initial unperturbed energy of 
the atom (see rescaling details below; atomic units are 
used) . It is also assumed that loq > 1 and that the initial 
orbital momentum I < (S/cj)^/"^. For eg < Ec the elec- 
tron energy E only performs small oscillations around 
its initial value and therefore ionization is impossible in 
the classical system. Above the chaos border instead, 
the electron's energy increases in a diffusive way with a 
diffusion rate per unit time given by 



(AEl 
At 
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(2) 



This diffusion leads to electron's ionization after a typi- 
cal diffusive time scale to ~ Eq/De, with Eo = — l/2nQ 
being the initial energy. Such classical diffusive ioniza- 
tion requires many microwave periods {tDio /2t: ^ 1) 
and quantum interference effects can suppress this dif- 
fusion leading to quantum localization of chaos [^,^. 
Such dynamical localization of chaos leads to a quan- 
tum probability distribution /at exponentially localized 



in the number of absorbed photons — {E — Eq)/u), 
namely Jn oc exp(— 2|A^0|/£0^). For the general case of 
monochromatic field excitation in a complex spectrum 
the localization length , measured in the number of 
photons, can be determined via the one-photon transi- 
tion rate F ~ Db/cj^ and the density of coupled states 
Pc §: 



2iTTpc = 3.33- 



,10/3 ■ 



(3) 



The last equality in the above equation corresponds to 



the hydrogen case, where pc 



— 'm3 



due to Coulomb de- 



generacy (which, as is known, is responsible for the ap- 
pearence of an additional integral on motion ||]). Quite 
obviously, in case of strong localization, namely when 
the localization length is less than the number of pho- 
tons required for ionization, TV/ = |i?o|/w = no/2wo, the 
quantum ionization is exponentially small and therefore 
negligible compared to its classical value. In the opposite 
case > TV/, namely above the quantum delocalization 
border e^o 



,7/6 
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(4) 



quantum ionization takes place in close to the classical 
case We note that the above dynamical localization 
represents a deterministic analog of the Anderson local- 
ization in disordered quasi-one-dimensional chains. In 
our case the site index corresponds to the photon num- 
ber Ntf,, while TV/ plays the role of the effective sample 
size. 

The above theoretical results have been checked by 
different groups and more recently reconfirmed in 
numerical simulations with newly developed algorithms 
The predictions of dynamical localization have 
been also confirmed by laboratory experiments [ni]-n3|. 
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However, in the mesoscopic solid-state language, the ef- 
fective "sample size" available in these experiments, cor- 
responding to Nj w 20 was too short to test with suffi- 
cient accuracy the theoretical predictions. In fact in this 
situation one observes significant mesoscopic fluctuations 
in the ionization border which have been discussed in de- 
tail in 1^,^ . While the main structure of these mesoscopic 
fluctuations can be well described by the quantum Kepler 
map it is still highly desirable to have much longer 
"samples" with Nj > 100 to study experimentally the 
dynamical localization of chaos in more detail. We note 
that recently the dynamical localization has been also ob- 
served in experiments with cold atoms propagation [p^ ; 
however, in this situation there are other experimental 
restrictions. 

In order to have larger Nj values, one has, either to 
increase hq or to take ujo 1. However, the first possi- 
bility is restricted by experimental conditions where one 
prefers to have uq < 100. The second choice leads to 
a regime close to the static field ionization in which no 
classical chaos exists; in any case, localization of classi- 
cally chaotic motion does not take place for << 1. 
In order to actually have large Nj values, it is neces- 
sary to take a different approach and work with atoms 
which are classically chaotic already in the absence of 
the microwave field. There are two main possibilities 
to have chaotic Rydberg atoms. One way is to put hy- 
drogen or Rydberg atoms in a magnetic field. In this 
case, the classical dynamics becomes chaotic when the 
Larmor frequency is comparable with the unperturbed 
Kepler frequency [p^-p^. Another possibility is to use 
Rydberg atoms with quantum defects in a static electric 
field. Recent investigations have shown that for a suf- 
ficiently strong static field the level spacing statistics is 
similar to the case of random matrix theory (RMT) [p^ . 
From the experimental viewpoint the case with a static 
electric field is simpler to deal with and in fact can be 
studied in laboratory experiments similar to [^-^ . For 
such atoms the chaos border for the microwave field drops 
to zero and therefore even at very small microwave field 
one can expect to see diffusive excitation in energy. In 
addition, such a diffusion can take place for much lower 
frequencies with l/ng <ti loq <ti 1, and this fact allows to 
increase the values of Nj up to few thousands. 

The investigation of the interaction of chaotic Rydberg 
atoms with a microwave radiation is also interesting from 
another viewpoint. Indeed, it is known that a chaotic 
structure of eigenstates leads to a strong enhancement 
of the interaction. In nuclear physics, as was shown by 
Sushkov and Flambaum [^, this effect leads to an en- 
hancement of weak interaction and parity violation by 
a factor of thousand or more. Also in the problem of 
Anderson localization in disordered solid-state systems 
it has been found that a short range repulsive/attractive 
interaction between two particles can strongly enhance 
their propagation p3|. All this indicates that a chaotic 



structure of Rydberg atoms can strongly increase their 
interaction with radiation. This should lead to a signif- 
icant decrease of the quantum delocalization border as 
compared to the standard case of internally non-chaotic 
atoms studied in [pUq]. In fact, the experiments by Gal- 
lagher et al. p^pl| indicated a lower ionization border 
than for the hydrogen case ||^. The physical mechanism 
of such ionization, proposed by Gallagher et al. pO| , ^ , 
is based on a picture of successive Landau-Zener cross- 
ings in a slowly oscillating microwave field. However, the 
question how such transitions can proceed to high levels 
was never studied in detail. Moreover, for ui > \/n^ one 
enters in a niultiphoton regime which was not discussed 
in (20|j2^. 

Due to the above reasons, the investigation of chaotic 
enhancement of microwave ionization of Rydberg atoms 
allows to address a new physical regime in which thou- 
sands of photons are required to ionize one atom. 
Our theoretical and numerical results indeed clearly de- 
mostrate the existence of such enhancement and provide 
a clear physical picture of the ionization process. In sec- 
tion II we discuss the case of atoms in parallel magnetic 
and microwave fields, while the case of Rydberg atoms 
in static and microwave fields is analysed in Section HI. 
The main results are discussed in Section IV. Some re- 
sults have been presented in pS-M. 



II. THE HYDROGEN ATOM IN MAGNETIC 
AND MICROWAVE FIELDS 



The Hamiltonian of a hydrogen atom in parallel mi- 
crowave and uniform magnetic fields writes 



2 



2p2 



1 



(5) 



+ ezsin(wi). 



where z is the direction of the fields, p — [x^ + j/^)^/^, 
lol = B / c = B{T)/Bq is the cyclotron frequency, Bq = 
2.35 X lO^T, e and uj are the microwave strength and 
frequency respectively (atomic units are used). Due to 
the cylindrical symmetry, the z component of the angu- 
lar momentum is a constant of the motion and here 
we consider — 0. 

The above Hamiltonian, expressed as a function of 
cartesian coordinates {xi}, their conjugate momenta 
{pi}, the time t and the parameters Lz, col, e, w has 
the property 



A-^e, \-^/^uj) = \-^H{{x,}, {p,}, t, Lz,ujL,e, w). 



(6) 
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If we choose A = l/rxg, with Eq = — l/2ng the initial en- 
ergy, the classical dynamics depends on no only via the 
scaled variables 



Lz/na, 



LUq 



(7) 



The coordinates and momenta scale as Xi = Xi/n^, 
Pi — PiTiQ and the rescaled time t — t/riQ, up to a factor of 
1 /27r, counts the number of Kepler periods in the electron 
motion in the absence of external fields (lul — e = 0). 

In order to study the dynamics of the time-dependent 
Hamiltonian (|^), it is possible to introduce an extended 
phase space in which the Hamiltonian becomes conserva- 
tive with respect to a fictitious time rj. The new Hamil- 
tonian writes K = H + ujo^j with (j) = ui^t = Lot and 
N as new classical conjugate variables. The equations of 
motion for 4> and N are given by 

d<j) _dK _ dN _ dK _ dH 

Notice that rj is equal to the scaled time t up to an ad- 
ditive constant, and therefore the ordinary Hamilton's 
equations follow for the other canonical variables. Since 
_ftr is a constant of motion, in the quantum case the 
change of N would give, apart from a l/ng scaling fac- 
tor, the number of photons AA^ exchanged by the atom 
with the field [|: AN = -AH/luq = -(l/no)(A£;/tj) = 
-(l/no)A7V. 

The singularity of the Hamiltonian (||) at r = (p^ + 
2;2)i/2 = can be removed by introducing the semi- 
parabolic coordinates u = (f -I- z)^/^, v — (f — z)^/^ and 
the regularized time a, defined as 



drj — dt ^ [u^ 



I da = 2rda, 



(9) 



which changes faster than i near the nucleus and more 
slowly far from it. The equations of motion generated by 
the Hamiltonian (||) are then equivalent to the equations 
of motion generated, with respect to the new time a, by 
the scaled and regularized Hamiltonian 



(10) 



Indeed, for any classical quantity f(u^v,(f);pu,Pv, N), we 
get 



df df dr] r,^ 2 , 2n 
da drj da 



(11) 



where {K, /} denotes the Poisson bracket between K and 
/. Thus, if we take Nq = — (1 / uq) {Eq / uj) as initial con- 
dition for the compensated energy K is equal to zero 
and we obtain 

{/C, /} = {u2 + v^J}K + {u^ + v^){K, /} = (12) 

that is /C works as an effective Hamiltonian for the time 
a. For = we have 



pI+pI ^ , (wLno)^2 2/ 2 , 2^ 



-2 



-U V [u -f f ) + 



(13) 



^ (u* _ sin (j) + u;oN{u^ + v^). 



For e = 0, uiqN —1/2 and therefore the Hamiltonian 
(|l3| ) represents two harmonic oscillators of unit frequency 
coupled by the term ^ {ujLnQ)'^u'^v'^ [v? -t-u^) which orig- 
inates from the diamagnetic interaction. For low scaled 
magnetic field w^ng, the motion is regular and the orbits 
are quasiperiodic. However, the diamagnetic term has 
cylindrical symmetry since it depends only on the per- 
pendicular distance p from the magnetic field axis. On 
the contrary, the Coulomb term has spherical symmetry. 
When these two terms are of comparable magnitude, 



1 



1 2 2 



(14) 



(that gives oJLn^ of the order of unity), then the com- 
petition between different symmetries leads to chaotic 
motion In Figj^ the classical phase space struc- 

ture is illustrated by the Poincare surfaces of section. For 
^L^o — 3 some islands of stability still exist but their size 
is small, while for lolu^ — 9.2 no regular structures are 
visible and the classical motion is dominated by global 
chaos [|6|-@. 

Due to this internal chaos, the turn on of the microwave 
field immediately leads to diffusive energy excitation of 
the electron, with classical diffusion rate per unit time 
Db =< {AE f > /At. From Eq.(|) E = ezuj cos(ujt) , 
with the typical frequencies for the z motion of the order 
of the Kepler frequency lok = (for ulUq not too large 
with respect to 1). For ^ 1, -Ds can be estimated in 
quasilinear approximation |p5| , giving Db ^ {ezoj)'^ /luk 
and then 



(15) 



where xi is a constant to be numerically determined and 
Dq = e^no/2 is the diffusion rate for ujl = and wq = 1 
when the microwave intensity is strong enough to induce 
chaos 

For Wo 3> 1, in analogy with the case without mag- 
netic field, the microwave interaction is mainly effec- 
tive when the electron passes close to the nucleus, where 
the Coulomb term dominates the diamagnetic term, and 
therefore, as for the case ujl = 0, one has P5[: 



Df 



X2-D0W0 



-4/3 



(16) 



with a constant X2 again to be numerically determined. 

An example of diffusive energy excitation, for scaled 
frequency ^ 1, is shown in Fig.|^. Here wltzq = 3, 
so that the motion is chaotic even in the absence of the 
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microwave and therefore a diffusive process occurs even 
if the microwave intensity is very smaU. 

The frequency dependence of the diffusion rate is 
shown in Fig.^, for wltt-q = 3 and 9.2. The asymptotic 
behaviors (both for smaU and large frequencies) are in 
agreement with the theoretical estimates ( p^ ) and ( p^ ) 
respectively, indicated in the figure by the straight lines, 
with the coefficients xi X2 numerically determined. 
These coefficients depend weakly on ujl, for magnetic 
field strong enough to induce internal chaos (wing ^ 2). 
Below this value the internal motion becomes integrable 
and for low enough eq diffusion drops to zero, as demon- 
strated in the insert of Fig. I (see also ||). 

The classical diffusive process will lead the electron to 
ionization in a time tj « Eq/Db, with Eq = — f/2nQ the 
initial electron energy. In the quantum case, for uq ^ 1 
the time evolution initially follows the classical diffusion. 
The excitation proceeds via a chain of one-photon tran- 
sitions, which eventually brings the electron into contin- 
uum. However, if the ionization time is sufficiently large, 
quantum interference effects may suppress the diffusive 
excitation leading to dynamical localization in the num- 
ber of photons. 

In order to check this possibility we numerically an- 
alyzed the quantum dynamics p6t, following the wave 
packet evolution in the eigenstates basis of the magnetic 
field problem with e = 0. 

To obtain these eigenstates at e = in a given energy 
window around the initially excited level with eigenen- 
ergy E\g w Eq, we diagonalized the Hamiltonian in a 
parabolic Sturmian basis |^,^. This basis is well suited 
since it is complete and discrete and can efficiently re- 
produce both the bound and continuum states of the hy- 
drogen atom. In addition, all the Hamiltonian matrix 
elements can be expressed in a simple analytical form 
and strong selection rules exist for the parabolic quan- 
tum numbers rti, n2 (An^ = 0, ±I,±2 for i = 1,2). A 
minor disadvantage of dealing with a Sturmian basis is 
associated with its nonorthogonality. Due to this fact, 
we had to solve a generalized eigenvalue problem of the 
type Hipi = Ex- Sipi , where also the overlap matrix S has 
strong selection rules (An^ = 0, ±1 for i = I, 2). As a 
result the matrices H and S are very sparse. To study 
the time evolution we used up to 1250 eigenstates of the 
above generalized eigenvalue problem. In order to obtain 
a good convergence of these 1250 eigenvectors and eigen- 
values we had to diagonalize matrices of size larger than 
13, 000. For these reasons the use of an efficient Lanczos 
algorithm has been crucial. The time evolution was 
computed by a split-step metod, similar to the one used 
in |]. 

In our computations, we considered as initial state an 
eigenstate at e = with eigenenergy E\g ^ Eq — — l/2nQ 
and the time evolution in the microwave field was fol- 
lowed up to r = 200 microwave periods. The pa- 
rameters were varied in the intervals 0.05 < ujq < 3, 



0.002 < eo < 0.04, 40 < no < 80, for ujLnl = 3 and 9.2. 

In this quasiclassical regime the quantum energy ex- 
citation has initially a diffusive character (see the insert 
of Figji]) and therefore it is possible to compute a quan- 
tum diffusion coefficient Dg^ from a linear fit, for the few 
first microwave periods, of the wave packet energy square 
variance < (AE)'^ > vs. time. Fig.^ shows that quan- 
tum and classical diffusion coefficients are similar over a 
wide range of frequencies (0.05 < ojq < 3 for uq = 60). 
This demonstrates that quantum dynamics mimics, for a 
finite interaction time, the classical excitation. 

In the quasiclassical regime, the one-photon transition 
rate Tb can be related to the classical diffusion rate |^ as 
Tb = Db/uj^- Indeed the change in energy produced by 
a one-photon transition is AE — ±a; and Tb measures 
the number of such transitions per unit time. Therefore, 
the ratio Tb/Tq, where Tq = Do'^o transition rate 

for the chaotic case at ujl = and = 1, is equal to 



Fb 
Fo 



DqUJ^ 



(17) 



This result is remarkable as it relates the quantum tran- 
sition rate to a classical characteristic of motion, namely 
the diffusion coefficient. In order to check the above esti- 
mate, Tb was numerically evaluated according to Fermi's 
Golden Rule: 



= ? E kmn|'{^(^^A,„ - Sa„ + ^)- 



S{Ex,„ - Ex„ - u;)}, 



(18) 



with z„i„ matrix elements of z between the eigenstates 
at e = with eigenenergies Ex^ , i?A„ • For the numeri- 



cal computation, the Dirac delta function in Eq.(18) was 
substituted by a Lorentzian function: S{x) — > f{x) = 
(l/7r)(e / (x^ -t- e^)), with e of the order of the level spac- 
ing. The validity of the classical-quantum correspon- 
dence ( |l7|) is confirmed in Fig.^, for uq — 40, 60, 80 and 
0.05 < Wo < 3. 

Even if initially the interaction with the microwave 
field results in a quantum diffusive excitation, quantum 
interference effects can suppress the diffusive behavior 
before ionization takes place. In this case, the diffusive 
broadening of the quantum probability distribution over 
unperturbed levels stops. The corresponding localization 
length (measured in number of absorbed photons) is 
proportional to the one-photon transition rate F b and to 
the density of coupled states pB- is = 2tiTbPb |§- For 
the chaotic case at lo^ — and luq = 1 the localization 
length is £^ — 3.3eo"-o ~ 27rZ?o"-oPo — SttFoPo, where 
Po — "-0 is the density of coupled states Actually, 
without the magnetic field there is an additional approx- 
imate integral of motion, related to Coulomb degeneracy, 
and therefore the density of coupled states pq is by a fac- 
tor no smaller than the number of levels per unit energy 
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interval nfj |^] . Therefore from the above expressions and 
Eq.0 we get MM 



Db Pb 



(19) 



This resuh provides a theoretical formula for the the 
quatum localization length, which involves only classical 
characteristics of motion, namely the classical diffusion 
coefficient and the density of coupled states. The latter 
is related, in the quasiclassical regime, to the phase space 
structure (see below Eqs.(|o|)-(|2^)). The estimate (|w|) is 
valid for £5 > 1, namely in the quasiclassical regime, in 
which a large number of photons is absorbed and a large 
number of levels is excited. Moreover, the microwave fre- 
quency should be larger than the average level spacing 
(ujpB > 1); otherwise, levels would move adiabatically in 
time leaving no room for diffusive energy growth . 

For wltt-q sufficiently large, the internal motion is 
chaotic: as a consequence, Coulomb degeneracy is re- 
moved and the density of coupled states pB ^ tiq- More 
precisely, pB = 0.34nQ (for lolu^ — 3, uq — 60) and 
PB = 0.14n4 (for uluI = 9.2, uq = 60) (see FigJ|). The 
density follows the dependence oc This is re- 

lated to the fact that the diamagnetic term is identical to 
a two-dimensional harmonic oscillator in the x — y plane. 
Indeed for a harmonic oscillator the spacing between lev- 
els is proportional to the frequency (here represented by 
lul) and therefore their density scales with clj]^^ ■ This re- 
sult is confirmed by a quasiclassical computation of the 
density of states (see the insert of Fig.|^). Here the num- 
ber fl(E) of quantum states having energy less than E is 
related to the corresponding volume in phase space 



n{E) 



1 



(2^)2 



dpdzdppdpz, 



H<E 



(20) 



where the Hamiltonian H is taken at e = 0. The density 
of states is 



PB 



dn{E) 4 dVl{E) 



(21) 



dE " dE 
where the volume in the scaled phase space is defined as 
1 



n{E) 



(27r)2 JH<B=-k 



dudvdpudpv, 



(22) 



with 



H 



pI + pI 



2(m^ + v'^) V? + 



u^v'^. (23) 



According to the results of Fig.^ such semiclassically de- 
termined density of states is in a good agreement with the 
results obtained by direct diagonalization of hydrogen in 
magnetic field. 

According to the estimate ( pj| ) and differently from the 
zero magnetic field case, the localization length is non 



homogeneous in the number of absorbed photons. In- 
deed, from equations (|l^) and ( |l6| ) we get Ib ^ n}^ ^ 
(7V/-7V^)-"/2 for Wo < 1 and^B - no - {Ni^N^)-^''^ 
for 3> 1 p5[ , with Nj = 71^/2^1^ the number of pho- 
tons required for ionization and = N — Nq, where 
A'o — Eq/uj. However, for £b <C Nj this effect is not im- 
portant and the quasistationary distribution is exponen- 
tially localized in the photon number N^. As for the case 
of hydrogen in a microwave field only , the probabil- 
ity distribution over unperturbed levels (e — 0) displays 
a chain of equidistant peaks in energy. The probability 
amplitudes in a one-photon interval decay esponentially 
with the photon number: ip^ oc exp(— |A'^ — Nq\/£b)- 

The value of the localization length is strongly en- 
hanced compared to the length = S.Segno/wQ*'^'^ j|] 
at wl — and ujq > 1. Actually, for a strong enough 
scaled magnetic field. Coulomb degeneracy is removed 
and the eigenfunctions ?/;, when developed on the basis 
of the hydrogenic eigenfunctions f, have a large num- 
ber of randomly fluctuating components, increased by a 
factor M = Pb/ Po' 



M 

fc=i 



(24) 



with Cik ~ 1/vM due to normalization condition. Be- 
cause of this, the external microwave perturbation V cou- 
ples the eigenstates with typical interaction matrix ele- 
ments 



V^nt = {'>pJ\V\^|^^ 



M 
k,l=l 



c*iCik{(pi\V\ipk) 



(25) 



Due to selection rules for V in the hydrogenic basis, Vint 
is the sum of order M uncorrelated terms and is therefore 
reduced by a factor 1 / ^/M in comparison with the case 
without magnetic field, whereas the spacing Ai? between 
nearest neighbors levels in the spectrum scales as 1/M. 
Thus, the admixture factor 



77 ~ V^nt/^E . 



'M . 



(26) 



is strongly enhanced, for high ng values, in comparison 
with the zero magnetic field case. As a result, the local- 
ization length £b , which is proportional to , is increased 
p6[ by a factor 



X2PbIP0 ^ M r^UQ. 



(27) 



A similar effect was discussed for two interacting parti- 
cles in a disordered solid state system : the interaction 
allows the two particles to propagate coherently on a dis- 
tance much larger than one-particle localization length. 
The statistical enhancement effect is quite general and 
takes place in different pysical problems: actually this 
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effect was first studied for the weak interaction parity- 
breaking in the scattering of polarized neutrons on com- 
plex nuclei Q . In the compound nucleus the level spac- 
ing is typically A_E « 1 eV, whereas the typical energy 
scale due to the strong interaction is AEq « 1 MeV. As a 
consequence, the number of principal components of the 
eigenf unctions in the one-particle basis is M « 10^ and 
therefore the admixture factor between states of opposite 
parity is enhanced by a factor ^/M « 10"^ compared to a 
usual parity breaking admixture of the order 10~^. 
The condition £b = Nj allows to determine the quan- 



tum delocalization border |25 2i 



no y 6.6Db Pb(^o 



(28) 



For eo > localization effects become unimportant and 
diffusive ionization close to the classical case takes place. 
For ti>o ^ 1 this border is (x2?^o/'^L«o)^^^ times smaller 
than the quantum delocalization border at zero magnetic 
field 1^ eqQ — uj^^ /^/6.6nQ. For wq ^ 1 the border is 
given by — (1^^71^/6.6x11^0)^ ^'^ /^o^ which remains well 
below the ionization border for a static electric field in 
the presence of a parallel magnetic field. Since the lo- 
calization length is non homogeneous in the number of 
absorbed photons and becomes larger near the ionization 
border, the actual value of the delocalization border 
will be slightly decreased as compared to (p8|). 

The above dynamical localization theory was tested 
by detailed numerical simulations of quantum evolution. 
The probability distribution f\ over the eigenstates at 
e = for the localized case is shown in Figs.|^,|| as 
a function of the number of absorbed photons = 
{E\ — Eq)/lu. In order to suppress fluctuactions this prob- 
ability was averaged over 20 microwave periods. To find 
the numerical value £b n of the localization length we first 
computed the probability /jv in each one-photon inter- 
val around integer values of N^. Then the least square 
fit with /at ~ exp(— 2iV0/^sAf) for > allows to 
determine £bn- The quantum distribution at luq — 0.2 
(Fig.^(b)) starts to deviate from the exponential profile 
for Nij, > 60 and a plateau appears. This can be related 
to the fact that localization length is non homogeneous 
in the number of absorbed photons {£b ^ (Nj — N^)~^^^^ 
for ujQ 1). In the figures also the classical distribution, 
normalized to one-photon interval, is shown. This distri- 
bution was obtained by numerical solution of Hamilton 
equations with 1000 trajectories initially distributed mi- 
crocanonically on the energy surface at energy Eq. The 
comparison of the probability distributions at different 
times (80 < T < 100 in Figs .0(a) J|(a) and 180 < r < 200 
in Figs. 0(b), ||(b)) shows that the quantum case is well 
localized over a quasistationary distribution. Its profile 
does not change significantly by increasing the interac- 
tion time, apart from a slight increase in the tail. We 
attribute the decrease in the probability decay at large 



Ncf, > to the growth of the diffusion rate Db when 
approaches 1 from below (see FigJ|). 

These distributions, together with the time depen- 
dence of the square variance of the photon number < 
(AN^)"^ > (see Fig.^) demonstrate how the quantum- 
mechanical energy diffusion is strongly suppressed in 
comparison with the classical case. 

The numerically obtained localization lengths £bn are 
presented Fig.|l^, for two fixed values of eo and different 
luq in the range 0.05 < < 3. Since no — 60, the number 
of photons required for ionization goes from Nj = 10 to 
Nj — 600. The straight line, corresponding to £bn = £b, 
demonstrates the fairly good agreement between numer- 
ical data and dynamical localization theory over a wide 
range of parameters (the localization length changes by 
two orders of magnitude). 

Finally, Figs. plJT^ show that the excitation at w^tiq — 
3 is much stronger than at zero magnetic field, due to 
the chaotic enhancement of electron's interaction with 
the microwave field. Since the condition for the transi- 
tion to chaotic dynamics at zero magnetic field is given 
by HJI eo > Cc = V49cj^^^, the case of Fig.|n| (that cor- 
responds to Fig. 10 in |^) is above the chaos threshold 
(eo > ec — 0.015). However the quantum distribution 
at = clearly demonstrates exponential localization 
of diffusive excitation since eo < e^o = 0.14. On the 
contrary, for lultiq = 3 we are above the delocalization 
border (eo > eq — 0.016) and thus numerical results show 
a good agreement between classical and quantum distri- 
butions. The quantum and classical square variance of 
the photon number are shown in Fig.|l^. In the absence 
of the magnetic field, the localization of the quantum 
motion after a few microwave periods is evident; on the 
contrary, for wltiq = 3 the quantum square variance fol- 
lows the classical one and starts to deviate only for r > 20 
because of the finite size of the basis. 



III. RYDBERG ATOMS IN STATIC ELECTRIC 
AND MICROWAVE FIELDS 

Another way to study the microwave interaction in an 
intrinsically chaotic atomic system is to consider alkali 
Rydberg atoms in a static electric field; this case is inter- 
esting as it is more suitable for laboratory experiments. 

The Hamiltonian of a hydrogen atom in a static elec- 
tric field is separable in parabolic coordinates and thus 
the classical motion is integrable. On the contrary, in 
alkali atoms the valence electron moves in an attractive 
potential in which the 1/r dependence is modified near 
the origin due to the charge distribution of the inner elec- 
trons, which form a core of size a few Bohr radius. Due 
to precession in a static electric field, the electron always 
collides with the core. This fact prevents the separation 
of variables and the motion becomes chaotic for certain 
field and energy values Q . 
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The potential within the many-electron core region 
can be reasonably approximated by a one-electron spher- 
ically symmetric model p9|: 



Vc{r) 



1 



exp(— air) — 012 exp(— asr) 



(29) 



where Z is the nuclear charge and the parameters 
(i = 1,2,3) are numerically determined so as to repro- 
duce quantum defects, that characterize the ionic core 
of alkali atoms in quantum mechanics. In this way, the 
overall potential V{r) — ~l/r + Vc{r) complies with the 
two conditions V{r) ^ —Z/r (for r — > 0, where the core 
fails to screen the nuclear charge) and V{r) ^ —1/r (for 
r — > (XI, where due to the screening the valence electron 
feels only a unit charge). Even if different core potential 
forms could reproduce quantum defects, both classical 
and quantum properties of the system are sensitive only 
to the size of these defects, not to details of the core. 

In the quantum case, the core potential Vc{r) intro- 
duces a phase shift ttJ; in the radial wave function. The 
electron energy is then given by Eni = —l/2{n — , 
namely the effective principal quantum number becomes 
n* = n — 61, where the quantum defect di depends on 
the orbital momentum / but changes only weakly with 
the energy E. Since the size of the core is small only 
the quantum defects for low angular momentum values 
{I < 3) are significantly different from zero. 

At the ionization threshold, the quantum defects Si are 
given, in quasiclassical approximation, by the action dif- 
ference in the nonhydrogenic and hydrogenic radial mo- 
tion between turning points ||30|| : 



Si — lim 

TT R~*oo 



2r2 



dr 



(30) 



''0 



1 (/ + 1/2)2 



2r2 



dr 



In this formula, l{l + 1) has been replaced by {I + 1/2)^ 
so that quasiclassical approximation remains valid also 
for small angular quantum numbers |]3l|] . The turning 
point for the hydrogenic part is tq = {I + 1/2)^/2, while 
for the nonhydrogenic part we obtain tqc as the root of 
Vi^Oc) + (l + l/2)2/2ro^ — 0. The parameters ai were 
fixed in order to reproduce the most important quantum 
defects So,Si,S2- For Li we assumed Sq = 0.4, Si — 0.04, 
S2 = 0; for Na So = 1.35, = 0.85, S2 = 0.01; for Rb 
So — 3.14, Si — 2.65, S2 = 1.35; we always considered 
Si>2 = [|2|. Typical values of the ai parameters are 
quite different from those given in Tab.l of j2^, probably 
because in that case these parameters were optimized so 
as to directly reproduce the energy levels. In any case 
this difference is not important in evaluating the quan- 
tum defects, that can be reproduced from Eq. (|30|) with 
ai parameters taken from p9] with a negligible error. 



The classical Hamiltonian for alkali atoms in parallel 
static electric field eg and microwave field esui{ujt), both 
in the ^;-direction, writes as 



p2 1 

H = ^ ^Vc{r)+esZ + ezsui{ijjt). (31) 

2 r 



The scaling property expressed by Eq.(^) remains 
valid, provided that also the static field and the core po- 
tential parameters are scaled as eso = Cs^oi "^io = otin^ 
(i = 1, . . . , 3). In other words, differently from the case 
of hydrogen, the problem depends separately on and 
7io and not only on the scaled parameter eso • As a matter 
of fact, the transformation ai — > aio = aiU^ changes the 
rescaled core size. 

The study of the classical dynamics of the above prob- 
lem was done in a similar way to the case of the hydro- 
gen atom in magnetic and microwave fields. The system 
was studied in an extended phase space that includes the 
time as a new canonical variable and the Coulomb sin- 
gularity at r = was removed by the introduction of 
the semi-parabolic coordinates u, v and the regularized 
time a (Eq.(^). For simplicity's sake, we considered only 
the case with z-component of the angular momentum 
L, =0. 

For e = 0, due to the core effects, we have a transition 
from regular to chaotic motion by increasing the static 
field tao ||l9|. However, the nature of chaotic dynam- 
ics in this case is qualitatively different from the case 
of diamagnetic hydrogen. In the latter case, the irreg- 
ular motion is due to the competition of the spherical 
Coulomb interaction, that dominates near the nucleus, 
and the cylindrical diamagnetic interaction, that domi- 
nates at large distances. In the former case instead, the 
chaotic motion arises from the presence of a core, which 
is small compared to the typical size of Rydberg atoms. 
Fig.^ shows that the phase space remains dominated by 
orbits which jump between different tori of the Coulomb 
problem. A trajectory follows a torus of the hydrogenic 
system until it encounters the core region where it is scat- 
tered on a different torus. As a result a single trajectory 
is able to explore most of the phase space energetically 
accessible to it. The lobes near v — Qin the Poincare sur- 
face of section in Fig.^ are due to the attractive core, 
which allows larger values of the momentum for a given 
energy. 

This chaotic motion affects also the quantum energy 
level spacing statistics ||3^. For integrable systems the 
level spacing statistics has generally the Poisson form 
P(s) — exp(— s), with s the spacing between two con- 
secutive levels, normalized to the local mean level spac- 
ing. For systems that display chaotic classical motion 
the distribution P{s) is described by the Random Ma- 
trix Theory (RMT). For a number of chaotic systems, 
with time reversal symmetry, the energy level statistics 
has been found in agreement with the predictions of the 
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Gaussian orthogonal ensemble of random matrices, with 
P{s) given by the Wigner-Dyson distribution 

P(s) = ^seM-js')- (32) 

A Wigner-Dyson distribution of level spacings has been 
seen in the Li atom in a static electric field strong enough 
to induce classical chaos . In Fig.|lJ the normalized 
nearest neighbor distribution P{s) is shown for esO = 0.02 
at uq = 60, in the range 55 < n < 72. The hydrogen dis- 
plays nearly Poissonian behavior while for Rb, in which 
quantum defects are large, P{s) is close to the RMT re- 
sults. The case of Na demonstrates a weaker level re- 
pulsion than Rb, while Li, which has smaller quantum 
defects, presents a distribution rather close to the Pois- 
son case. 

The energy levels for alkali atoms in a static electric 
field (the so-called Stark levels) were obtained by matrix 
diagonalization in a spherical hydrogenic basis; quantum 
defects were introduced simply in the alkali energy levels 
(for angular momentum Z < 3), while dipole matrix ele- 
ments were taken as in the hydrogen case. This approach 
is valid outside the ionic core, where V{r) is essentially 
hydrogenic and so is particularly well suited for Rydberg 
states, in which the electron mainly remains far from the 
nucleus. The off-diagonal dipole matrix elements for the 
z operator decrease rapidly with the difference in energy 
between initial and final states. As a result. Stark eigen- 
values and eigenvectors in a given energy window can 
be obtained by diagonalization of Hamiltonian matrices 
of limited size around this region. A limitation of the 
method is that only bound states are included in the ba- 
sis set and increasingly large matrices must be used to 
obtain convergence as the static field ionization border 
Es'T-^ = 0.13 is approached. 

The number of mixed hydrogenic states can be char- 
acterized by the inverse participation ratio 

^x. ^ iJ2\^xAnJtr\ (33) 

n,l 

where {n, I) are the coefficients of the expansion of the 
Stark eigenstate with energy E\. on the spherical basis. 
Fig.^ shows how the inverse participation ratio ^ (ob- 
tained averaging in one-shell intervals E{n — 1/2) < 
< E{n + 1/2), with E{x) = -l/2x^) increases with 
the energy for a given static field or equivalently (in the 
insert) with the field for a given energy. As a result, more 
and more shells are mixed moving to the saddle point 
Eso = 0.13 and chaotic properties become stronger. For 
Eso = 0.02, no = 60, ^ w 75 and this indicates that eigen- 
functions are significantly spread over different shells. 

When the microwave field is turned on, the internal 
chaos leads to diffusive energy excitation in the classi- 
cal dynamics, with energy diffusion rate per unit time 
Dc =< {AEy > /At. As for the hydrogen atom in 



static electric and microwave fields, the rate Dc can be 
compared with the diffusion rate Dq for hydrogen with 
only the microwave field present and at loq — 1. Again 
the asymptotic behaviors Dc = Xi-Do'^o (^^r loq <^ 1) 

and Dc = X2Dq'^o (^°^ ^ 1) are expected. This 
is confirmed by Fig.|l^, that shows the frequency depen- 
dence of the scaled diffusion rate Dc/Dq for Rb, Na and 
Li at eso = 0.02, eo = 0.005, hq = 60. The core poten- 
tial parameters ai are not changed with the initial level 
no- As a result, the dependence on the initial energy 
or on no can no longer be eliminated by classical rescal- 
ing. However the classical scaling rule remains a useful 
approximation since Dc/Dq varies very slowly with no 
for Rydberg atoms (for Rb at e^o = 0.02, eo = 0.005, 
luq — 0.1, Dc/Dq remains between 0.019 and 0.023 for 
20 < n < 80). This fact can be understood by taking 
into account that only for low I values the electron col- 
lides with the core. The energy diffusion rate is therefore 
determined by the frequency of these collisions, which is 
determined by the frequency of classical precession in /, 
namely by the scaled Stark frequency Wso — 3eso, inde- 
pendently from no. In other words, the dimensions of 
the core play no important role for the classical diffusion 
rate. 

Fig.^ presents some other interesting features: a 
sharp resonant peak for luq = 1 and a plateau for 0.02 < 
ujQ < 0.4. This latter characteristic is due to the fact that 
core collisions, that are responsible for energy diffusion, 
occur at the Stark frequency Usq, independently from 
ujQ. This fact is also illustrated in the insert of Figjl^: 
for very low microwave intensity (eo — 0.0005 <C e^o) 
and for ujq = 0.1 the diffusion coefficient grows with the 
static field eso, while for eo = 0.005, when eso < coi the 
microwave field dominates the angular momentum pre- 
cession and the dependence of Dc/Dq on e^o becomes 
rather weak. 

The peak at wo = 1 in Fig.|l^ is due to the resonance 
between Kepler and microwave frequency. Figjl^ shows 
that the diffusion rate drops abruptly as the electron goes 
out of the resonance (for eo = 0.005 this happens after ap- 
proximately 10 microwave periods). Such sharp change 
of the diffusion rate was observed only near the main res- 
onance luq = 1, while for other integer values of ojq our 
studies did not indicate any sharp change. 

To study diffusion in the quantum case we consid- 
ered the initial wavefunction at t = as an eigen- 
state of the Hamiltonian ( |3T| ) at e = with eigenenergy 
E\o Eq — — l/2ng As in the case of the hydro- 

gen atom in magnetic and microwave fields, the quantum 
dynamics was followed by a split-step method, similar to 
the one described in |^ . The quantum diffusion rate Dq 
was obtained by a linear fit of the energy square variance 
< (AE)'^ > vs. time for the first few microwave periods. 
The classical-quantum comparison for the scaled diffu- 
sion coefficients Dc/Dq and Dq/DQ (see FigJT^) shows a 
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quantitative difference. This is due to the fact that only 
low angular momentum values contribute to the diffusion 
and so purely quantum effects can be expected. In addi- 
tion, the influence of the ionic core is stronger in classical 
mechanics because trajectories can penetrate in arbitrar- 
ily small regions in the phase space. On the contrary, 
quantum mechanics tends to smooth over such regions. 
Apparently this is the reason for which classical diffu- 
sion rate is sistematically larger than the corresponding 
quantum value. 

Due to quantum interference effects, the above dif- 
fusive excitation may eventually stop before ionization. 
The resulting quasistationary distribution is character- 
ized by an exponential decay in the number of absorbed 
photons, with a localization length £g proportional to the 
one-photon transition rate Tq and to the density of cou- 
pled states pc- iq = i-KVqPc [||. The rate is given by 



and since the internal motion is chaotic the 
density of coupled states is pc = n^- This leads to |27|: 

D„ 



^1 ^^-PTrh^o, 



(34) 



where, as in Eq.(|l9|), = 3.3eQnQ and the conditions 
£q > 1, uipc > 1 must be satisfied. 

In order to check the theoretical prediction (|3^), the 
quantum evolution of an eigenstate with eigenenergy 
E\g — l/2riQ, and no — 60 was followed up to 200 
microwave periods in the Stark basis for eso = 0.02. The 
total basis size was up to 1150 states. The system pa- 
rameters were varied in the intervals 0.003 < Cq < 0.03, 
0.02 < Wq < 2, which corresponds to 15 < Nj < 1500. 

Typical examples of stationary distributions for Rb at 
luq ~ 0.08 and no — 60 are shown in Fig.|l^. In or- 
der to determine the numerical value IqN of the localiza- 
tion length, we first computed the total probability /at 
in each one-photon interval. Then the least square fit 
with with ''^ exp(— 2iV0/^gAr) allows to determine £qN. 
In Fig.|l^ the numerical localization length iqjy agrees 
with the theoretical £q value (Eq.(|3|)) for eq = 0.005. 
The figure also shows how quantum excitation is strongly 
suppressed in comparison with the classical case. For a 
smaller microwave intensity {cq = 0.002) the fit in the 
interval < < 20 gives igN = 3.7, whereas the fit 
for the tail {N^ > 50) shows a much slower decay, with 
iqN — 28. We checked that this slope change is not af- 
fected by the variation of the basis size and integration 
step. This change in the probability decay at large A^^ 
can be attributed to a significant change in the eigenstate 
structure in a static electric field when one approaches 
the saddle point behind which tunneling takes place. In- 
deed according to the results shown in Fig.|l^ a Stark 
eigenstates for higly excited levels projects on a larger 
number of hydrogenic levels. Additional investigations 
are required for a better understanding of this effect. 

The comparison of the numerically obtained localiza- 
tion lengths £qN vs. the theoretical estimates iq con- 



firms the predictions of Eq.(|3J) (see Fig.^ for Rb and 
Na at a fixed eo value, with uiq varied over the wide range 
0.02 < uJo < 0.5). We did not study dynamical locahza- 
tion for the Li atom since it is close to the integrable 
case. 

Equation (pM|) allows to determine the quantum dclo- 



calization border |27| from the condition £q = Nj: 



6.6Da 



(35) 



For £0 > diffusive ionization takes place. Actu- 
ally, this threshold should be lowered since ionization is 
also possible due to the static field above the threshold 
A more accurate estimate is given by the 



relation £„ 



N^, with Aff = Ni{l 



Also one 



should keep in mind that the estimate ( |35| ) is based on 
the initial local value of Dq/Do taken at rt « tiq. In the 
process of excitation the ratio Dq/Dg can be changed, for 
example due to a sharp peak near ~ 1 (see Fig.p^. 
This can give some additional decrease for the border e^. 

Fig.|2l| shows that above the delocalization border 
(eo > eq) the quantum excitation is close to the clas- 
sical one. In this case the wave packet escapes into con- 
tinuum and quantum interference effects are unable to 
freeze quantum diffusion before ionization. We note that 
the border (^5|) is lower than for the hydrogen atom p7| ] 
approximately by a factor ^/nQ which appears due to in- 
ternal chaos originated by quantum defects (see also Fig. 5 
in (27)). 

The border (|3^) for alkali Rydgerg atoms in static and 
microwave fields is in qualitative agreement with a series 
of experiments by Gallagher et al. ||2^,^, that showed 
at low frequency (wo <C 1) a scaled microwave ionization 
threshold eg ~ I/tiq instead of the static field hydrogenic 
border e^o — 0.13. Also for Rb atoms the ec threshold is 
well below the one for hydrogen jl^. Using the data for 
Dq in the frequency interval 0.02 < u)q < 0.5 one can see 
that the quantum delocalization border is approximaely 
tq w 1.5eG (see Fig.6 in ||). 

The ionization process was interpreted by Gallagher et 
al. [ po|j2l| as due to a chain of Landau-Zener transitions 
to higher-lying states, until the static field ionization bor- 
der is reached. Indeed, in the presence of quantum de- 
fects, the first avoided crossing between the n and n + 1 
Stark manifolds occurs at a field eo = l/3no ~ to- How- 
ever, this theory doesn't explain which is the dynamical 
mechanism that brings the electron through the whole 
chain of Stark levels up to ionization. In addition, the 
Landau-Zener theory does not apply if the microwave 
frequency is larger than the typical spacing between lev- 
els, namely for wo^-o > 1- On the contrary, the dynam- 
ical localization theory allows to understand ionization 
of atoms in a static electric field for the non-adiabatic 
regime ojo^-o > 1- According to this theory the mecha- 
nism of ionization is qualitatively different from the one 
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proposed by Gallagher et al, namely ionization takes 
place due to diffusive excitation in energy originated by 
internal chaos existing in absence of the microwave field. 

A difficulty for the direct comparison with experiments 
[ p^ , pO 21 3^ is that only few of them were done in the 
presence of a a static electric field and, in addition, 
they were in the regime ujpc — loqUq < 1. However, there 
is a case for Na at hq — 28, initial orbital momentum 
la = 2, e^o = 0.024, loq = 0.027 (Fig.2d in |2|]) which is 
not far from our conditions {ioorio — 0.76). The exper- 
imental 50% ionization threshold is eoex — 0.002. This 
value is about 20 times smaller than the quantum delo- 
calization border given by Eq.(^5|) with Dq/Do = 0.0027. 
Indeed, quantum simulations in the eigenstate basis, ex- 
tended up to 200 periods (see Fig.|2|), give a localization 
length £qN = 4 <C iV/ = 520. We understand this dis- 
crepancy as due to deviations in the tail of probability 
distribution {ig^ — 25 for Jn < 10^^) similarly to the 
case discussed in Fig.|l^. 

In order to support this interpretation, we followed the 
time evolution in the hydrogenic basis and simulated tun- 
neling ionization by an absorption mechanism. Namely, 
we modified the evolution of the wave function ip{t) as 

i;{t) = Teyi-p{-iHt) exp(-|i)V'(0), (36) 

where T is time ordering operator, 7 is a diagonal op- 
erator in the spherical basis, with matrix elements 7„ w 
l/2TTn^ for all the states in a shell with principal quan- 
tum number n > (es^-s = 0.13). In this way, there 
is absorption for levels above the static field ionization 
threshold after a time approximately equal to the Kepler 
period T„, = 27rrig. In this model, ionization probability 
grows with time in a nearly linear way (see Fig.|23|), with 
ionization rate (per microwave period) T = 2.5 x 10^*. 
This fact can become important for very lon g in teraction 
times (r/ = 8.3 x 10^ microwave periods in ||2l|]) leading 
to strong ionization. 

Another situation in which the experimental ioniza- 
tion border is much lower than the one given by the 
theoretical estimate ( |35| ) was observed in Beterov et al. 
experiments |]35| . In this experiments, with no = 36, 
Iq = 1, Wo = 0.51, e^o = 0, approximately half of Na 
atoms were ionized after t/ = 1.5 x 10^ microwave peri- 
ods at Co — 0.003. A static field of strength ego = 0.02 
would give a diffusion rate Dq — 0.057I?o- In such a case 
the quantum delocalization border is eq = 0.032 which is 
much higher than the experimental value. The numeri- 
cal simulations with effective absorption discussed above 
give an absorption rate F = 9.3 x 10"'^ for eq = 0.003 
and r < 200. Such ionization rate would give a signifi- 
cant ionization during the long interaction time tj used 
in experiments | ]35[ |. 

A possibility, alternative to the presence of a weak 
static electric field, is that some noise existing in the 
waveguide could destroy localization and give a larger 



ionization probability compared to the theoretical expla- 
nation (^). 

For cjQ <C 1, one could expect that the slowly varying 
microwave field will also play the role of a static electric 
field even if — 0. However, in reality, this expectation 
can be valid only if luq is much smaller than the frequency 
ujso of classical precession determined by the Stark split- 
ting and equal to 3eo {^a ^ 3eo). If this condition is not 
satisfied then the mixing of I states will not take place 
and this is in agreement with our numerical data for Rb 
at no = 60, Iq = 0, wq = 0.1, eo = 0.01, Cso = 0, where 
after r = 100 microwave periods only few Z-states are 
mixed (< ^ 2.8). At the same time in a presence of 
a static field e^o = 0.02 the probability spreads over all 
accessible Z-values (< I >= 28). The localization in I 
space for Rb atoms was found also in |Q at Cg = 0. 

The existing experiments does not allow unfortunately 
to make a direct check of dynamical localization theory 
because the interaction times were too long and therefore 
it is difficult to control the effect of environment. How- 
ever present laboratory conditions allow to study short 
interaction times (r sa 200) and high quantum numbers 
no ~ 60. In these conditions, according to our numerical 
and theoretical results, quantum excitation of atoms is 
well described by the dynamical localization theory. Ex- 
periments in this regime will allow to test the quantum 
localization effects in a range of parameters much larger 
than it was so far possible. 

Additional investigations should be done for the regime 
ujQ > 1. In this case the electron precess with very low 
frequency tOso ^ ^^o a-nd therefore the static electric field 
gives an adiabatic perturbation on the localized distri- 
bution in the photon number. The case of such adia- 
batic destruction of the localized case was discussed in 
1^ and manifestations of this effect for alkali Rydberg 
atoms should be analyzed more carefully. The situa- 
tion for alkali Rydberg atoms in a static electric field 
is quite different from the case of the hydrogen atom in 
magnetic and microwave fields, where for o^o > 1 the 
microwave frequency is of the order of the Larmor fre- 
quency {lulo ~ '-^o)- On the contrary, for alkali atoms 
LUso ^ ^^0 E^nd therefore localization takes place faster 
than the spreading of the wave function over the whole 
energy surface. 



IV. CONCLUSIONS 

In this paper we have analyzed the properties of mi- 
crowave ionization of chaotic Rydberg atoms. Similar 
to the cases of parity violation in nuclei |^ and of two 
interacting particles effect in disordered systems |2^], a 
chaotic structure of eigenstates (in absence of microwave) 
leads to a chaotic enhancement of radiation interaction 
with atoms. As a result the localization length in the 
number of photons is strongly increased as compared to 
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the usual situation in which the internal dynamics of the 
atom, without the microwave field, is integrable ; as 
a consequence, the quantum delocalization border drops 
down significantly. The theory of dynamical localization 
developed for such chaotic atoms is in good agreement 
with the results of extensive numerical simulations. 

Investigations of such atoms in laboratory experiments 
represent a new important opportunity to provide de- 
tailed results for quantum chaos and dynamical localiza- 
tion. Indeed, due to internal chaos, the excitation pro- 
ceeds in a diffusive way even if the microwave frequency 
is much smaller than the Kepler frequency of electron's 
rotation (wq ^1). As a result the number of photons Nj 
required for ionization can be as large as few thousands, 
thus allowing to investigate the dynamical localization of 
chaos in great detail. Experimental conditions are rather 
similar to those in previous experiments |]l3| , |l9| - p2|j35| 
and are available now in modern laboratories. 

Here we have discussed the case of parallel fields in 
which the magnetic quantum number Lz remains an inte- 
gral of motion. It can be also interesting to study a more 
general situation with arbitrary field's orientation and 
polarization. In this case the density of coupled states 
will be even larger (pc ~ ng). This can lead to an addi- 
tional growth of localization length by a factor uq and to 
a decrease of the delocalization border by a factor ^/no- 
At the same time the effective "sample" size Nj can be 
increased by ng for such small frequencies as utQ ~ l/f^o• 
However, this case deserves more detailed studies. In- 
deed, the existence of additional approximate integrals 
of motion or the appearance of very slow adiabatic fre- 
quencies is not excluded (especially for a static electric 
field) and this may lead to new interesting results. 

The results of the present paper also show that theoret- 
ical and experimental studies of chaotic Rydberg atoms 
still represent a challenge for fundamental research of 
quantum chaos. 
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FIG. 1. Poincare surfaces of section at scaled energy 
E = H{e = 0) = -1/2 and (a) uluI = 3, (b) uluI = 9.2. 
The sections are on the v—pv planes defined by u = 0, Pu > 0. 
In (a) five regular and one chaotic trajectories are shown, in 
(b) only one chaotic trajectory is followed. The centre of the 
atom corresponds to i; = 0. 




FIG. 2. Example of diffusive energy excitation 

{{AEf)/E^ with the number of microwave periods r, for 
uiltio = 3, eo = 0.005 and ujo ~ 0.1. The straight line 
fit allows to determine Db/Do = 0.062 via the relation 
< (A£)2 {2-KDB/uj)r. 




UJQ 

FIG. 3. Dependence of the scaled diffusion rate Db/Do 
on the scaled frequency ujo for eo = 0.005, aJLTig = 3 (full cir- 
cles) and LOhTi^ = 9.2 (open circles). The st raig ht lines show 
the theoretical dependence for ci;o ^ 1 (Eq.(|l5[)) and ci;o 3> 1 
(Eq.(|l|)), with XI = 25, X2 = 0.8 {uLul = 3), xi = 16, 
X2 = 1.6 [ujLnl = 9.2). Ensembles from 200 to 1000 trajecto- 
ries, initially distributed microcanonically on the energy shell, 
have been used. The insert shows Db /Do as a function of the 
scaled Larmor frequency WiWo, for eo ~ 0.005 and ujq = 0.1. 
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FIG. 4. Classical (full circles) and quantum (open trian- 
gles) scaled diffusion rates, for LOLriQ — 3. In the quantum 
case no = 60, 0.01 < eo < 0.04. The insert shows an ex- 
ample of diffusive energy excitation in the quantum case, for 
i^LTio = 3, eo = 0.02, Luo — 0.1; the straight line fit gives 
Db/Do = 0.055. 



FIG. 6. Scaled density of states ps/n^ as a function of 
n for LdLnf-i — 3 (full circles) and uJLnQ — 9.2 (open circles) 
at no — 60. Solid and dashed lines refer to the quasiclas- 
sical computation (Eq.(^l|)). The insert gives quasiclassical 
scaled density of states vs. scaled Larmor frequency. The 
straight line shows the theoretical asymptotical dependence 
(ps/no — c/ujLriQ, with c — 1.5). 




FIG. 5. Gamma ratio Fb/Fo vs. scaled frequency ujo, for 
lolti^ = 3, no = 40 (full line), no = 60 (dotted line), no — 80 
(dashed line). The one-photon transition rate Fs is evalu- 
ated according to Fermi's Golden Rule (Eq.(p^), with the 
Dirac delta function substituted by a Lorentzian function, 
with spread e — 0.5/ ps, where pB = 0.34no is the density 
of states. Full circles show classical values of Db/Dqujq at 
UIlUq = 3. 
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FIG. 7. Probability distribution vs. number of absorbed 
photons N^: quantum distribution fx over the eigenbasis at 
e = (full curve); quantum probability in one-photon inter- 
vals /jv (circles); classical probability in one-photon intervals 
(dotted line). The straight line shows the fit for the expo- 
nential localization. Parameter values: no = 50, = 1-3, 
^LiT-o — 3, eo ~ 0.008; scaled diffusion rate Db/Dq — 0.53, 
theoretical localization length estimate (Eq.(^9|)) £b = 2.8; 
(a) 80 < r < 100, £bn = 2.1; (b) 180 < r < 200, £bn = 2.4. 
In the classical simulation 1000 trajectories are included in 
the microcanonical ensemble. 
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FIG. 8. Same as in Fig. 7, with no = 50, uo — 0.2, 
LULul = 3, eo = 0.007, Db/Dq = 0.090, £b = 15.4, (a) 
80 < r < 100, £bn = 18; (b) 180 < r < 200, £bn = 18.8 
(in both cases the fit is limited to < < 60). The classi- 
cal ensemble consists of 1000 trajectories. 
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FIG. 9. Dependence of the quantum (full curve) and 
classical (dotted curve) photon number square variance 
< (AA'^^)^ > on the number of microwave periods r for 
no — 60, LOo — 0.1, eo — 0.005, loltiq = 3. The classical 
ensemble consists of 5000 trajectories. 



FIG. 11. Chaotic enhancement of microwave excitation. 
Same as in Fig.^ with no — 66, ~ 2.5, a^^no — 3, 
eo = 0.04 > = 0.016, Db/Dq = 0.97, 40 < r < 50. FuU 
triangles, linked by a dashed line, give quantum probability in 
one-photon intervals at the same parameter values but with- 
out magnetic field, 290 < r < 300 and with initial value of 
the orbital momentum quantum number Iq — 5. In the latter 
case, Ec = 0.015 < eo < e^o ~ 0.14, theoretical — 1.1, 
numerical —1.7. 
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FIG. 10. The numerically computed localization length 
£bn vs. the theoretical estimate £b (Eq.(p^)) for eo = 0.005 
(full circles) and eo — 0.01 (open circles). Data are for 
LJLng = 3 or 9.2, no = 60, 0.05 < loq < 3. The straight 
line corresponds to £bn = £b ■ 
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FIG. 12. Quantum (full curve) and classical (dotted 
curve) square variance of the photon number vs. the num- 
ber of microwave periods for the case of Fig.^ Above and 
upper scale: loltiq = 3, below and lower scale: u>Lno — 0. 
The case lol = is followed for a longer interaction time than 
for loltiq = 3 to emphasize dynamical localization effects. 
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FIG. 13. Poincare surface of section for Na at eso ~ 0.02, 
no = 60. The v — pv plane of section is defined by u = 0, 
Pu > 0. The section is constructed following only one trajec- 
tory. Core potential parameters are ai — 2.48, 02 ~ 0.54, 
as = 1.43. 



FIG. 15. Inverse participation ratio ^ (Eq.(p3|)) for Rb, at 
EsO ~ 0.02 and no — 60, as a function of n. The insert shows 
the dependence of ^ on eso, for no = 60. 
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FIG. 14. Level spacing statistics P{s) for alkali atoms in 
a static electric field, with eso = 0.02, no = 60, 55 < n < 72, 
for (a) H, (b) Li, (c) Na and (d) Rb. The dotted line shows 
the Poisson distribution, the solid line is the Wigner-Dyson 
distribution (Eq.(||)). 
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FIG. 16. Dependence of the scaled classical diffusion rate 
Da/ Dq on the scaled frequency uiq for alkali Rydberg atoms 
in a static electric field, at eo = 0.005, tao = 0.02, no = 60 
for Rb (circles) , Na (triangles) and Li (squares) . The straight 
lines show the theoretical dependence for ujq and ojq ^ 1, 
with xi = 140, X2 = 0.3. Ensembles from 200 to 1000 tra- 
jectories, initially distributed microcanonically on the energy 
shell, have been used. Core potential parameters for Na as 
in Fig|l|, for Rb ai = 3.10, ^2 = 1.56, Q3 = 0.76 and for 
Li Qi = 2.84, Q2 = 4.05, 0:3 = 4.26. The insert shows, for 
Rb at uja = 0.1, no = 60, the dependence of Dc/Dq on £^0 
for eo = 0.0005 (diamonds) and the more flat behavior for 
eo = 0.005 (circles). 
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FIG. 17. Dependence of ((A£)^) /El on the number of mi- 
crowave periods r, for no — 60, eso = 0.02, eo = 0.005, uo = 1- 
The straight line fits give Dc/Dq = 1.4 (for < r < 10) and 
D^/Dq = 0.17 (for 10 < r < 70). 
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FIG. 18. Quantum scaled diffusion rate Dq/Do vs. fre- 
quency tJo for Rb (full circles) , Na (full triangles) and Li (full 
squares) at 6^0 = 0.02, no = 60, 0.005 < eo < 0.03. Open 
circles show the classical scaled diffusion rate Dc/Do for Rb, 
with parameters oii as in Fig.[L6l 
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FIG. 19. Probability distribution vs. number of absorbed 
photons A^^: quantum distribution fx over the Stark eigen- 
states (full curve), quantum probability in one-photon inter- 
vals /iv (circles) and classical probability in one-photon in- 
tervals (dotted line) for Rb at no = 60, tao = 0.02, ujo — 0.08, 
eo = 0.005, 180 < r < 200. Here DJDo = 0.0047 and 
the theoretical localization length (Eq.(p4)) is £q — 13. The 
straight line shows the fit for the exponential localization, 
with £qN = 13. Diamonds give /jv (shifted down by 10^) for 
Rb under the same conditions, except for eo — 0.002, iq = 2.1. 
The straight hue shows the fit for < iV,;, < 20 {£qN = 3.7), 
whereas the fit for A^^ > 50 gives £qM ~ 28. Classical data 
are obtained with parameters as in Fig. 16. 
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FIG. 20. The numerically computed localization length 
£qN VS. the theoretical estimate £q (Eq.(|3^)) for Rb (full cir- 
cles) and Na (full triangles) at no — 60, eso = 0.02, eo = 0.005, 
0.02 < u)o < 0.5. The straight line corresponds to IqN = £q- 
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FIG. 21. Same as in Fig.^ but above the delocaliza- 
tion threshold tq (Eq.(|5|)), for Rb at no = 60, e^o = 0.02, 
LJo = 0.08, DJDo = 0.0047, eo = 0.03 > = 0.027, 
40 < r < 50. 
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FIG. 22. Same as in Figjl^, for the conditions of Fig. 2d in 
[21], for Na at no = 28, £^0 = 0.024, ujo = 0.027, eo = 0.002, 
180 < r < 200, Dg/Do = 0.0027, £q = 1.1, fqiv = 4 (straight 
line) for /jv > 10"^ and E^n = 25 for /at < 10^^ 
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FIG. 23. Probability distribution over the spherical hy- 
drogenic basis for the case of Fig.^. Time evolution is fol- 
lowed under the same static and microwave field conditions, 
starting from a hydrogenic eigenstate with no = 28 and initial 
value of the orbital momentum quantum number lo = 2, with 
absorption parameters (Eq.(^)) given by 7„ = 1/2-Kng for 
Esn* > 0.13 (n > Us = 42). The insert shows that ionization 
probability depends on time in a nearly linear way (for the 
interaction time here considered, r < 200), with ionization 
rate (per microwave period) F — 2.5 x 10^*. 
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